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ABSTRACT 

We have carried out a high-sensitivity search for circumstellar disks around Herbig Be stars in the 
continuum at 1.4mm and 2.7mm using the IRAM interferometer at the Plateau de Bure (PdBI) . In 
this letter, we report data on three well studied BO stars, MWC 1080, MWC 137 and R Mon. The two 
latter have also been observed in the continuum at 0.7 cm and 1.3 cm using the NRAO Very Large 
Array (VLA) . We report the detection of circumstellar disks around MWC 1080 and R Mon with 
masses of M^ ~ 0.003 and 0.01 Mq, respectively, while for MWC 137 we estimate a disk mass upper 
limit of 0.007 Mq. Our results show that the ratio M^/M* is at least an order of magnitude lower in 
Herbig Be stars than in Herbig Ae and T Tauri stars. 

Subject headings: Radio continuum: stars - Circumstellar matter - Stars: individual (MWC 1080, 
MWC 137, R Mon) - stars:pre-main sequence 



1. INTRODUCTION 

The existence of accretion disks around massive stars 
(M> 5 Mq) remains a matter of debate. There is in- 
creasing evidence for the existence of flattened structures 
(disks) around high-mass protostars. However there is no 
clear evidence of disks in late r phas es, namely in Herbig 
Be (HBe) stars. iNatta et all 1)20001) compiled the inter- 
ferometric observations in mm continuum around HBe 
stars and found that the occurrence of disks in HBe is ~ 
0. This contrasts with the case of the lower mass HAe 
stars (2Mo < M, < 5 M©) in which the occurrence of cir- 
cumstellar disks is similar to that in T Tauri stars. They 
propose that the lack of disks around HBe stars is due to 
the rapid evolution of these objects, whic h disperse the 
surro unding dust and gas in about 10^ yrs ijFuente et al.l 
120021) . However, this lack of detection can be a sensi- 
tivity effect. HBe stars are usually further away than 
HAe stars, and higher sensitivity interferometric obser- 
vations are required to detect circumstellar disks around 
these objects. In this Letter, we present the first results 
of a high-sensitivity search for circumstellar disks around 
three BO stars, MWC 1080, MWC 137 and R Mon, which 
are the best studied of the five Herbig BO stars in the 
northern sky hsted bv lThe etlHI l|l994|) . 

2. OBSERVATIONS 

Interferometric observations in the continuum at 
1.4mm and 2.7mm have been carried out towards the 
Herbig Be stars MWC 1080, MWC 137 and R Mon, us- 
ing the IRAM^ array at Plateau de Bure, France, in the 
CD set of configurations. MWC 137 and R Mon, were 
also observed with the NRAO^ Very Large Array (VLA) 
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at 0.7 cm and 1.3 cm in its D configuration. Flux cal- 
ibration is accurate within 10% at 2.7mm and 20% at 
1.4mm in the PdBI images and within 20% in the VLA 
images. No correction for primary beam attenuation has 
been applied. PdBI and VLA images are shown in Fig 1. 
Flux densities at the star positions are shown in Table 1. 

3. MWC 1080 

Although, we have not observed MWC 1080 at cm 
wavelengths, some informatio n abou t the cm emission 
can be found in the literature. Curiel et al. (1989), from 
10" resolution observations at 6cm, detected a weak ex- 
tended component which includes the star and the reflec- 
tion nebula. jSkinncr et al. (199j) with higher angular 
resolution (~ 1") tentatively detected at 3 .6cm a com- 
pact s ource at the star position. Recently, IGirart et alJ 
()2002l) detected three 6cm sources using the VLA with 
an angular resolution of ~ 5". One of these sources is 
coincident with the star position. They did not detect 
a counterpart of this emission at 2.0cm with an upper 
limit of < 0.32 mJy. 

Intense emission is detected in the millimeter contin- 
uum images observed with the PdBI (see Fig lab). In- 
terferometric observati ons at 2.7mm in MWC 10 80 were 
previously reported bv lDi Francesco et al.l l|1997j) with a 
non-detection. Our image improves by a factor of 10 
the sensitivity and by a factor of 2 the angular resolu- 
tion of their data. The millimeter continuum emission in 
MWC 1080 is extended. Two intense components con- 
nected by a bridge of weak emission are observed in the 
2.7mm image. These components are much closer to the 
star than the VLA 6cm sources detected bv IGirart et alJ 
(2002). We will refer to these components as "NW" and 
"SE" . The better angular resolution of the 1.4mm image 
resolve the "NW" and "SE" clumps in several compo- 
nents, "NWr',"SEl", "SE2" and "Star". The total flux 
in the 1.4mm image is ~ 24 mJy. In our single-dish data 
we measured a 1.4mm flux of 73.8 mJy at the star pos i- 
tion w ith an angular resolution of ^ 10'.' 5 |Puente_et_alJ 
119981) . Thus the interferometer has resolved the 1.4mm 
emission and missed ~ 70% of the flux. The " Star" com- 
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Table 1 
Observed fluxes. 
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Fig. 1. — Continuum images obtained with the PdB and 
VLA arrays. The star marks the optical position [MWC 1080 
: 23:17:25.574 60:50:43.34; MWC 137: 06:18:45.504 +15:16:52.4; 
R Mon: 06:39:09.94 +08:44:10.0] and the filled triangle the po- 
sition of the radio source (Skinner et al. 1993). Contour levels 
are: a. -2.0 to -0.2 mjy/beam and 0.2 to 4.0 mjy/beam by 0.2 
mjy/beam; b. -1.6, -0.7, 0.7 to 3.4 mjy/beam by 0.9 mjy/beam; 

0.2 to 2.2 mjy/beam by 0.4 mjy/beam; 
0.3 to 3.6 mjy/beam by 0.6 mjy/beam; 
-0.6 mjy/beam 0.6 to 6.6 mjy/beam by 
-0.2 mjy/beam, 0.2 to 2.4 mjy/beam by 
-0.5 mjy/beam, 0.5 to 4.0 mjy/beam by 
-0.2 mjy/beam, 0.2 to 0.8 mjy/beam by 
-0.2 mJy/beam 0.2 to 1.4 mjy/beam by 0.2 



c. -0.2 mjy/beam, 

d. -0.3 mjy/beam, 

e. -1.8 mjy/beam, 
1.2 mjy/beam; f. 
0.2 mjy/beam; g. 
0.5 mjy/beam; h. 
0.2 mjy/beam; 
mjy/beam. 



ponent is coincident with the optical star position and 
remains unresolved in the 1.4mm image. This implies 
that its size is < 2" which corresponds to a spatial extent 
of <2000 UA at the distance of MWC 1080 (d=1000pc). 
This size is compatible with the emission arising in a 
circumstellar disk. 

4. MWC 137 

MWC 137 is associated with the one-arcminute size 
HII region S266 ( Fich 1993). Observation s of this region 
at 3.6cm were reported bv ISkinner et all l(l993r ). Our 
1.3cm image of MWC 137 consists of a compact source 
coincident with the HBe star surrounded by a thin shell 
with a radius of 30" -40" from the point source. Weak 
emission is found towards the thin shell and the point 
source. In the 0.7cm image, the compact source appear 
more intense while the emission of the thin shell becomes 
fainter. Since our aim is the detection of a circumstellar 
disk, we are only interested in the compact component 
(see Table 1). 

The continuum millimeter images show intense emis- 
sion centered at the star position. The total flux in the 
1.4mm image is ~ 10 mJy, i.e. 30% of the flux niea- 
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sured in the single-dish observations ijFuente et al.ll 99^ . 
From our 1.4mm image we derive a source size of 1.8" 
X 0.8", which corresponds to « 2300 AU x 1000 AU 
assuming a distance of 1300 pc, and a position angle of 
75±15deg. Note that the accuracy of the position angle 
estimate is heavily limited by the S/N ratio of the data. 

5. R MON 

R Mon was observed at 6cm, 3.6cm and 2cm using the 
VLA bv Skinne r et al. (1993). They detected continuum 
emission at a position slightly shifted from the optical 
one. A compact source is detected in R Mon at 1.3cm 
and 0.7cm. The position of this source is offset (0.74",- 
0.45") from the optical position but coincident with the 
radio source. A jet-like feature is observed in the 1.3cm 
image which is surely dri ven by the radio sou rce. This 
jet was already detected bv lBrueel et alJ l)1984|) in optical 
lines. Since the uncertainty of the optical position is not 
known, the identification of the radio source with R Mon 
is uncertain. 

We have detected continuum emission towards the 
star in the 2.7mm image with a peak emission of 4.1 
mJy/beam. The total integrated flux is 6.4 mJy re- 
vealing that the emission is extended. Deconvolving the 
2.7mm image with the beam, we estimate an emission 
size of 3" -4", whic h correspond s to « ^ 3000 AU at the 
distance of 800 pc. iNatta et al.l 1)2000^') reported a total 
2.7mm flux of 13. Oil. 3 mJy towards R Mon based on un- 
published data by Mannings (1998). This flux exceeds by 
a factor of 2 our present result. This reinforces our con- 
clusion that the emission is extended. Then, the different 
beams and calibration uncertainity can account for this 
discrepancy. The 2.7mm continuum emission is elongated 
in the direction perpendicular to the bipolar nebula axis 
and the centroid of the emission lies to the northeast of 
the position of the radio-source. The shape and position 
of the 2.7mm emission suggest that the 2.7mm emission is 
tracing a different emission component from that traced 
by the emission at cm wavelengths, and argues in favor 
of the existence of a flattened structure (disk) around the 
star. 

6. SPECTRAL ENERGY DISTRIBUTIONS (SEDS) 

The SEDs in the cm-mm range of MWC 1080, 
MWC 137 and R Mon are shown in Fig. 2. In or- 
der to avoid the contribution of the envelopes of these 
sources to the continuum emission, only interferometric 
measurements have been included in these SEDs. Even 
though, some problems still remain. The emission at 
cm wavelengths in MWC 1080 is extended. The neg- 
ative index derived by comparing the flux at 6cm and 
3.6cm is very likely due to the larger beam of the 6cm 
flux ( 5.5" X 5.03") compared with that of the 3.6cm 
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Fig. 2. — Spectral energy distributions (SEDs) at cm and mm wavelengths of MWC 1080, MWC 137 and R Mon. The fluxes at cm 
wavelengths have been taken from Skinner et al. (1993) and Girart et al. (2002), and at mm wavelengths are from this paper. The 
short-dash lines are model predictions for the free-free emission component. The long-dash lines are model predictions for circumstellar 
disks of 0.003 Mq and /3=1 (MWC 1080), 0.007 Mq and /3=1 (MWC 137), and 0.01 Mg and /3=0.5 (R Mon). The solid line is the sum 
of both components. The short-long straight line in central panel is the predicted fluxes for the emission arising in a stellar wind with a= 
+0.76 in MWC 137. 



one (1.24"x0.83"). The disk emission should be un- 
resolved at the distance of our source. Therefore, we 
are not interested in the extended component and con- 
sider the 6cm flux as an upper limit. Unfortunately, the 
free-free emission arising in the stellar wind is also ex- 
pected to be unresolved with our beam and could make 
an important contribution to the continuum flux even 
at millimeter wavelengths. In order to derive the disk 
mass in MWC 1080, we have assumed that the 3.6cm 
flux is due to free-free wind emission and subtracted it 
from the observed mm fluxes. These corrected values 
are then interpreted as optically thin disk emission, ac- 
cording to the expression F\ — d~^MdH\B\{Td) with 
K — 0.01(A(mTO)/1.3)~^cm^'^g^^. We have assumed a 
constant value of Td = 215 K following the recommen- 
dation of Nat ta et all f200G) for BO stars. With these 
assumptions, we obtain from the 1.4mm flux a disk mass 
of 0.003 M0. Since we have only an upper limit to the 
2.7mm flux, the value of (3 cannot be determined. We 
have compared the disk mass derived with this simpli- 
fied expression to those of more elaborated models (see 
DuUemond et al. 2002), varying the disk outer radius 
and the surface density profile over a large r ange of pos- 
sible v alues, and confirmed the conclusion of lNatta et al.l 
l)2000(l that the values of derived in this way are ac- 
curate within a factor 2-3. We have also considered the 
possibility that the 3.6cm emission arise in an optically 
thin HII region instead of free-free emission. With only 
two points in the SED, we cannot discern between these 
two possibilities. In this case, we obtain a slightly larger 
value of the disk mass of, ~ 0.004 Mq. Therefore we 
estimate that the mass of the disk around MWC 1080 is 
0.003±0.001 Mq. 

The complete SED of MWC 137 can be fitted with a 
single component with spectral index a — -|-0.76±0.01. 
This spectral index is consistent with that expected in 
the free-free emission arising in the stellar wind. Al- 
though a value of a =+0.6 is expected for a ionized 
isotropic wind, small deviations of this value can be ex- 
plained by a different geometry or a partially ionized 



wind. However, we cannot discard the existence of a 
very low-mass disk. In order to estimate an upper limit 
to the disk mass, we have extrapolated the 3.6cm flux to 
millimeter wavelengths with an spectral index of +0.6. 
Then, we have estimated the disk mass from the excess 
flux at 1.4mm. With the same values of k and T^ than 
in the case of MWC 1080, we obtained an upper limit to 
the disk mass of 0.007 Mq. 

The SED of R Mon is the most complete one. We have 
been able to fit all the centimeter data using a spectral 
index, acm ~ 0.64, which is characteristic of a spherical 
stellar wind. In order to estimate the mass of the cir- 
cumstellar disk, we have extrapolated the free-free emis- 
sion to millimeter wavelengths and fit the millimeter part 
of the SED with a simple disk model. As in the case 
of MWC 1080, we have assumed an unresolved disk at 
a constant temperature, T^ = 215 K and adjusted the 
value of Mrf and (3. Only values of /? = 0.25 - 0.5 fit our 
observational points. A value of (3 larger than 0.5 would 
produce an excess of flux at 1.4mm. This low value of (3 
could be interpretred as an optically thick compact disk, 
or an optically thin disk with large grains. The first pos- 
sibility would imply millimiter fluxes larger than those 
measured. A small disk of R = 100 AU should have a 
mass of 0.35 M© in order to have Ti.4m„i~l. But in this 
case, the 1.4mm flux would be at least an order of mag- 
nitude larger than our upper limit. The situation would 
be worse for a larger disk. The second possibility seems 
more plausible. Similarly low values of (3 have been found 
in Herbig Ae and T Tauri disks and reveal highly pro- 
cessed dust (iBeckwith fc Sargentil99ltlTesti et air2 003). 
With the fiducial value of the dust emissivity, we derive 
a disk mass of 0.01 Mq. Note, however, that, if the low 
value of (3 is interpreted as evidence of grain growth, one 
needs to reduce t he fiducial 2.7mm dust emissivity by a 
factor of about 4 l|Testi et al.ll2003D . and to increase the 
estimated disk mass by a similar amount. 

7. DISK DETECTIONS IN HBE STARS 

Our results prove the existence of circumstellar dust 
around MWC 1080 and R Mon. However, it is very dif- 
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Fig. 3. — Plot of the and M^/M* versus the spectral type 
of the star (top panels) and versus the stellar age (bottom 
panel) for a wide sample of HAEBE stars. Data have been taken 
from the compilation by Natta et al. (2000) (empty symbols) and 
the results reported by .Fuente et al.. 1,2001) and this paper (filled 
symbols). Dots mean detections while triangles mean upper limits. 
Note that disk masses in HBe stars are at least an order of magni- 
tude lower than in HAe stars. 



ficult to interpret it beyond doubt as evidence of a cir- 
cumstellar disk. In general, a compact and small 1.4mm 
source (< 1" - 2") at the star position is considered as 
a disk detection. This criterium seems sufScient in TTS 
located at distances of ~ 150 pc, in which an angular 
resolution of 1" corresponds to a linear scale of 100 - 200 
AU , typical size of circumstellar disks. Since HBe stars 
are usually located at distances larger than 1 Kpc, even 
the higher angular resolution provided by interferome- 
ters, ~ 1", corresponds to sizes of a few thousands of 
AU, and it is not possible to distinguish between "bona- 
fide" circumstellar disks and a flattened structure of a 
few 1000 AU. 

Another difficulty comes from the fact that many HBe 
stars are likely associated to lower mass companions that 
will not be resolved in our beam. MWC 1080 is a clos e 
binary with a separation of 0.7 " l | Le inert et alJll997f) . 
Mid-infrared observations bv iPolomskT^r^lJ l)2002() us- 
ing the Keck II telescope resolve the binary and found 
that warm material (~ 350 K), with a 3 - 10 fim spec- 
tral index consistent with the existence of a circumstellar 
disk, is surrounding MWC 1080 A. A cooler dust compo- 
nent is associated with MWC 1080 B. We have detected 
an unresolved mm continuum source in MWC 1080. We 
consider that this source is a good candidate for disk de- 
tection around a HBe star. However, since the angular 
resolution of our observations is not sufficient to resolve 



the two binary components, we cannot exclude the pos- 
sibility that the emission is associated with the younger 
companion. 

R Mon is also a close binary with a separation of 0.69" . 
Its companion is a very young TTS. IClose et a,].l I^I997^ 
found an extinction of Av = 13.1 mag towards the star 
which they interpreted as due to an optically thick disk 
of R = 100 AU. Our 2.7mm image reveals the existence 
of a compact source of « 3000 AU x 2300 AU. This 
size is larger than that expected for a circumstellar disk 
and suggests the existence of a flattened envelope sur- 
rounding the disk. Our value of the disk mass may be 
overestimated due to the contribution of the flattened 
envelope. However, the analysis of the mm spectral in- 
dex shows that the grains responsible of the emission are 
heavily processed, so that most of the emission is very 
likely arising in a circumstellar disk. If this is true, it 
is unlikely that the disk could be associated to the TTS 
companion, since in this case one sho uld use in deriving 
the disk mass a value of Td - 15 K l|Natta et alJl200fy 
and the estimated mass would increase to values in the 
range 0.2-0.5 Mq, too large for a TTS disk. 

We have completed the compilation of (interferomet- 
ric) disk observat ions of HAeBe stars carried out by 
iNatta et all lj2000|) with the observations in MWC 1080, 
MWC 137 and R Mon reported in this paper, and the 
observations i n HP 200775 and LkHa 234 reported by 
iFuente et alJ 1)200 ID . In Fig. 3, we plot the disk masses 
as a function of the spectral type of the star and the stel- 
lar age. We have plotted MWC 1080 and R Mon as disk 
detections; however, as discussed in the previous section, 
it is possible that the "disk" emission is either overesti- 
mated or not associated to the Herbig Be star itself, but 
to unresolved companions, thus reinforcing our conclu- 
sions. In any case, disk masses in HBe stars are at least 
an order of magnitude lower than in HAe stars. Plotting 
Md/M* as a function of the spectral type, one finds that 
while the value of M^/M* is roughly constant and equal 
to 0.04 for stars with spectral type A0-M7, M^/M* < 
0.001 in HBe stars. The same effect remains when one 
plots the disk mass against the stellar age. The values of 
disk masses in HBe stars are systematically lower than 
in HAe stars for stars with ages between 10^ and 10^ yrs. 
This lack of massive disks in HBe stars even for stars as 
young as 10^ years put several constraints to the possi- 
bility of forming planetary systems around massive stars. 
In a time-scale of < 10^ years the disk mass around a HBe 
star is similar to or even lower than the total mass of the 
planets in the solar system 0.0013 Mq). Thus planet 
formation should occur very fast, for planetary systems 
to exist at all around massive stars. 
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